Europhysics Letters 



PREPRINT 



Mean first-passage times for an ac-driven magnetic mo- 
ment of a nanoparticle 

S. I. Denisov 1 ' 2 , K. Sakmann 1 , P. Talkner 1 and P. Hanggi 1 

1 Institut fur Physik, Universitdt Augsburg - Universitdtsstrafie 1, D-86135 Augsburg, 
Germany 

2 Sumy State University - 2 Rimsky-Korsakov Street, 40007 Sumy, Ukraine 



PACS. 05.10.Gg - Stochastic analysis methods. 

PACS. 75.60. Jk - Magnetization reversal mechanisms. 

PACS. 75.50.Tt - Fine-particle systems; nanocrystalline materials. 



Abstract. - The two-dimensional backward Fokker-Planck equation is used to calculate 
the mean first-passage times (MFPTs) of the magnetic moment of a nanoparticle driven by 
a rotating magnetic field. It is shown that a magnetic field that is rapidly rotating in the 
plane perpendicular to the easy axis of the nanoparticle governs the MFPTs just in the same 
way as a static magnetic field that is applied along the easy axis. Within this framework, the 
features of the magnetic relaxation and net magnetization of systems composed of ferromagnetic 
nanoparticles arising from the action of the rotating field are revealed. 



Introduction. - The mean first-passage time (MFPT), i.e., the average time that elapses 
until a stochastic process reaches a prescribed domain, is an important characteristic of the 
considered process. It is widely used for describing various dynamic features such as exit 
problems, activation rates, lifetimes of metastable states and other noise induced phenomena 
[1,2]. The class of stochastic processes for which the MFPT can be calculated explicitly is 
rather limited [1]. In fact, the most general analytical results were obtained for continuous 
one- dimensional Markov processes within an approach based on the backward Fokker-Planck 
equation [3-7]. It is important to note that most of the exactly known results are restricted 
to Markov processes that are homogeneous in time. 

However, Markov processes that describe time- dependent systems usually cannot be ap- 
proximated by homogeneous processes. Well-known examples are, the phenomenon of Stochas- 
tic Resonance [8-10] and directed transport in ratchet-like systems [11] where external time- 
dependent fields play a crucial role. Since the above-mentioned method is no longer ap- 
plicable, the development of new approaches for calculating MFPTs in periodically driven 
systems presents an important challenge. For slowly varying external driving forces adiabatic 
and semiadiabatic approximations are available [12,13]. The path integral formulation of the 
conditional probability provides a convenient basis for approximations if the noise is weak [14]. 

In this Letter we present an analytical approach to the two-dimensional MFPT problem for 
a rapidly driven magnetic moment of a ferromagnetic nanoparticle. It is based on the backward 
Fokker-Planck equation in a rotating coordinate system, which describes the homogeneous 
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dynamics of the magnetic moment, and on an averaging procedure [15], which is similar to 
that used for describing the Kapitsa pendulum with oscillating point of attachment [16]. The 
analytical results favorably compare with numerical simulations of the equivalent Langevin 
equation. 

Basic equations. - In spherical coordinates, the Landau-Lifshitz equation [17] for the 
magnetic moment m(t) = m(sin 9 cost/?, sin 9 sin ip, cos 9) of a single-domain ferromagnetic 
nanoparticle assumes the form 

f ^0^ + 4z)w, <P= — ly-fsinfl^-Al-W, (1) 



m sin 9 \ 89 dip J ' m sin 2 9 \ 39 dip 

where 9 = 9{t) and ip — <p(t) are the polar and azimuthal angles of m(t), respectively, 7(> 0) 
is the gyromagnetic ratio, A(> 0) is a dimensionless damping parameter, m = |m(t)| denotes 
the conserved total magnetic moment, and W — W(9,ip,t) is the magnetic energy of the 
nanoparticle. Besides the damping, the interaction with a heat bath also generates stochastic 
elements in the dynamics of the magnetic moment. These effects lead to a random contribution 
to the effective magnetic field which is assumed to be given by isotropic, Gaussian and white 
noise [18, 19]. The conditional probability density P(8, ip, t\8' , ip', t') of the resulting Markovian 
process satisfies the forward and the backward Fokker-Planck equations which are equivalent 
to each other [20]. At equal times P(9, (p, t'\8' ', ip', t') = 5(8 — 9')S(ip — p'), where S(x) denotes 
the Dirac 6 function. Here wc will use the backward equation, which propagates t' from t on 
backward in time. In the present case it takes the form [21] 



A 7 2 df 

where A = Afc^T '/7m denotes the thermal noise intensity, fee the Boltzmann constant, T the 
temperature of the heat bath, L + (t') the backward Fokker-Planck operator, 



and the functions 



^^ t ) = -A^( A8in ^ + 4 )lr 



(3) 



g(9, <p, t) = 1 , ( sin 9^- -X^-)W 
A7771 sin 9 V 08 dip 

are proportional to the drift components of the angles 9 and ip, respectively, see cq. ([1]). Note 
that if the energy W depends on time, the backward operator is also time dependent. 

For the calculation of first-passage times a domain Q(t) of the magnetization state space, 
which is a sphere with radius |m|, has to be specified. In general this domain may change 
upon evolving time. Starting out in f2(t') the magnetization eventually will leave this domain. 
The instant tf of first crossing the boundary dtt(t) defines the first-passage time tj — t' . In 
order to prevent the trajectory to recross this boundary, absorbing boundary conditions must 
be imposed at dQ(t). For the backward equation they read 

P{9,<p,t\9',ip',t') = for (9',ip') £dSl(t'). (4) 

The probability Pn(t\9' , <p',t') that the magnetization has been in the prescribed domain for 
all times, starting out at t' up to time t, i.e., m(s) £ fl(s) for all t' < s < t, follows from the 
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conditional probability in the standard way: 

P n (t\9',ip',t')= [ d9d<pP(9,<p,t\9',tp',t'). (5) 
Jn(t) 

The integral of P n (t\9' , tp', t') over all times t determines the MFPT T n = (t f - t'), 

POO 

T u (9',<p',t') = / dtP n {t\6',v',t'). (6) 



Using all, the definition ©, the condition that Pn(t'\9', ip', t') = 1, and the backward Fokker- 
Planck equation ([2]), we find that the MFPT satisfies the following equation: 

A 7 2 L+(t>)TU9',cp',t>) + -^T n (9',ip',t>) = -1. (7) 

Accordingly, one obtains as boundary condition from eq. (U]) 

T n {6',<p',t')=0 for (#V) e dtotf). (8) 

For a periodic time dependence (with period T) of the backward operator and the fi-domain, 
i.e., for L + (t + T) = L + (t) and Q(t + T) = £l(t), the asymptotic solution of eq. is also 
periodic in time with the same period T. For time homogeneous processes and domains it is 
constant with respect to time and fulfills the well known Pontryagin equation A^ 2 L + Tq = — 1. 

General results. - We now consider the magnetic moment of a nanoparticle with uni- 
axial anisotropy and assume that a static magnetic field H is applied along the easy axis of 
magnetization which we choose as z axis, i.e., H = (0,0, H). Additionally, a rotating mag- 
netic field h(t) acts perpendicular to this axis, i.e., h(t) = h(cosojt, ps'muit, 0), where ui is the 
angular frequency and p = — 1,+1 corresponds to clockwise and counterclockwise rotation, 
respectively. Hence, the magnetic energy of the nanoparticle is given by 

W = \mH a sin 2 9 — mH cos 9 — rah sin 9 cos tj> (9) 

with H a denoting the anisotropy field and ip = ip — put. We note that in this case a detailed 
analysis of the deterministic dynamics of m(i) is presented in [22]. According to eq. ([9]), the 
functions / and g in eq. ([3]) depend on ip' and t 1 only through the single variable ip' = ip' — put': 

f =-2a(cos9'+H) sin9'+ — (\cos9' cos tp'- sin ?/>'), 

2ah (10) 
g = —(cos6' +H)- - a — (Xsmib' + cos6''cosV''). 
A A sin 9' 

where a = mH a /2kBT is the anisotropy barrier height in units of thermal energy fcgT, 
H = H/H a , and h = h/H a . By introducing a rotating frame, in which h(t) = /i(l,0,0) and the 
azimuthal angle ip is replaced by ip, the time derivative d/dt' goes over into d/dt' — pud /dip', 
and consequently the backward equation for the MFPT reads 



d 2 T n 
d9' 2 



+ ^I'W 2 - + (cot 9+f) ^ + {9 - paUu}) W = r { + ~dF) ' ( u) 



where t r = 2/(XyH a ) is the characteristic relaxation time of the precessional motion of the 
magnetic moment. Moreover, we assume that the domain Q(t) is stationary in the rotating 
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frame, i.e., f2(i) is bounded by a curve 4>n(ip) on the sphere. Then, the steady-state solution 
of pip is independent of time and can be found as the solution of the stationary Pontryagin 
equation 

d 2 T n , 1 d 2 T n ^.^n,, . ,dT n 

-de^ + ^W^ + {cotd + f) ~dW + {9 ~ paUuj) W = ~ atr - (12) 

This presents already a major simplification of the original problem. However, still a partial 
differential equation in two spatial dimensions remains to be solved. 

In the following we will consider two types of domains which we distinguish by the index 
fl := ±1. The domain with Q = +1 (denoted as up domain) contains the up magnetization, 
9 = 0, and is bounded in the rotating frame by a curve <j)+i(ip). Accordingly, the Q = 
— 1 domain (down domain) contains the down magnetization, 9 — w, and is bounded by 
a curve (p-\(ip). It is convenient to present the respective MFPTs out of these domains 
in the form Tq = Tq(9') + Sq(0',i/j'), where the overbar denotes an average over ip', i.e. 
(•) = (l/27r) J Q 27r dip'(-). Substituting the expansions f = f + fi and g = g + gi into eq. (fT2)l . 
one obtains an equation for Tn, reading 



(™t9' + f)—+f 1 —+g 1 — = -at r (13) 



d9' 2 v J ' d9' • 89' * 

and an equation for Sn, reading 



PSn, 1 d 2 S ( , ( ^.^dSa dT a .dS Q dS n dS Q 

^ + ^^ +{coi6+f) ^ +h ^ +{9 - paUu) w~ h ^^ gi w^°- (14) 

We emphasize that these equations are exact, i.e., they follow from the backward Fokker- 
Planck equation ([3|) for stationary domains in the rotating frame. 

High-frequency limit. - In the case of an arbitrary rotating field we are not able to solve 
eq. (fj"2|) analytically. But in the asymptotic limit of a fast rotating field, i.e. for uj ^> uj r = ^H a , 
eqs. and lfT4"|) can be solved readily. The key observation leading to their solution is that 
Sn —> for u! — v oo. Assuming also that the derivatives of Sq tend to zero as 1/lu for lo — > oo, 
eq. (|14|) simplifies to read in this high-frequency limit 

dSn dTn n , x 

patrUJ W fl ^o T= ( 5) 

With /i = (2ah/ X)(Xcos9' cos-0' — sin?/;'), see eq. (fTO|). the solution of this equation that 
satisfies the condition Sn = can be written in the form 

Sn =P^ (Acos 9' sin + cos V')^. (16) 

where <D = uj/io r . Note, that according to the above assumptions, this solution is valid 
if h/Cj <C 1. The u> dependence assumed in the derivation of eq. (fl~5| is now confirmed 
selfconsistently. Using eq. (|16[) and the relation gi = — (2ah /X)(X sin ip' + cos9' cosip')/ sin9' , 
we proceed to calculate the averages 



(17) 
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with h c ff — — ph 2 /Co. Finally, substituting eq. (fTT)) and / = — 2a(cos 9' + H) sin 9' into eq. (|13p. 
we obtain the desired equation for Tq in the high-frequency limit: 

d 2 T dT 

— H + [cot^ - 2a(cos6>' + H + h eS ) sin 9'}—£ = -at r . (18) 
dt) z dt)' 

This equation exhibits the remarkable result that a magnetic field that is rapidly rotating 
in the plane perpendicular to the easy axis of the nanoparticle acts on the nanoparticle's 
magnetic moment precisely as a static magnetic field h c g (in units of H a ) which is applied 
along the easy axis. The direction of the effective field h c g and the direction of the field 
rotation follow the left-hand rule, and the value of h g is the same for up and down domains. 
It is important to emphasize that although the condition h/Q <C 1 holds, the effective field 
can be large if h 3> 1. 

We note that 9' = 0, ir are singular points of eq. (fTS)) for Q = +1, —1, respectively. At these 
points the general solution of eq. (fT5|) exhibits logarithmic singularities. To prevent this non- 
physical behavior of Tq the regularity condition dTn/d6'\gi= 7r (i-n)/2 = must hold [23]. In 
order to derive the boundary condition for eq. (f!8|) . we decompose the function 0n(VO m t° its 
average and its periodic parts, i.e., <f>n(ip) = 0q + ^n(V0- Assuming that |i?nWOI ^ 0Hj fr° m 
the absorbing boundary condition of the full two-dimensional problem, Tq(4'q(' 1 I ) ')>' i P') = 0; 
we find the periodic part of </>n(V0> ^n(V') = — ph(\ cos (f> n simp + cosip)/ui, and the desired 
boundary condition for the averaged MFPT, Tn((j) n ) = 0. The solution of eq. (fl8|) with the 
specified regularity and boundary conditions becomes 

_ r cos9 ' „-a(x+H+h c!! ) 2 r n _ _ 

Tn{d')=at r _ dx e / dye<y+ H+h ^ , (19) 

J cos (pQ i- X J x 

where 9' G [0, if £1 = +1, 9' e [</>_i,7r] if Q. = —1, and the angles 4>q can be chosen 
depending on physical situation. 

In the case of a high potential barrier, a ^> 1, and moderately large total effective fields, 
\H + h e e\ < 1, the magnetic moment resides near one of two equilibrium directions, up or 
down. Since transition times between these states by far exceed the relaxation times towards 
these states, the averaged MFPT Tn(9') describing the transition from one state f2 to the 
opposite state —Q only weakly depends on the precise location of the initial magnetization, 
as long as 9' lies within the domain of attraction of the considered state Vt. Also the precise 
location of the absorbing boundary 4> n has practically no influence on Tq(9') if it is located 
well beyond the separatrix which divides the state space into domains of attraction of the up 
and down magnetization. Then the eq. Q19[) yields in leading order in a 

fZ a[l+n(H+h ca )] 2 
T n =t r J = = = = . (20) 

V a 2 [1 - (H + /i cff ) 2 ][l + n(H + h cB )} 

Using eqs. (|16|) . (fl9|) and (f20|) . it is not too difficult to demonstrate that |Sh| <C |ATq|, where 
ATji = Tq — T 1 ji|^ =0 is the contribution of h to To- This means that the periodic part of 
Tn(6', ip') can be neglected such that Tn(9', ip') « Tq. 

In Fig. [T] the theoretical prediction (JTHJ) and the asymptotic approximation (f2T)]) are com- 
pared with the results of a numerical simulation of the coupled Langevin equations for the 
two angles 9 and cj> which are equivalent to the process described by the eq. ([2]). 

Next we examine the most interesting case of zero static field, H = 0, when only the 
effective field influences the MFPTs. Specifically, if \h eS \ < 1 then eq. $2® yields T n = 
tq exp(f22a/i c ff), where tq — (t r / 4) ^/ 'ir / a exp a. Thus, the rotating magnetic field increases the 
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Fig. 1 - (Color online) The natural logarithm of the dimensionless MFPT, T — T+i • ui r , is displayed 
as a function of the dimensionless anisotropy barrier height a. The solid line is the result of present 
theory (|19|) . and the broken line depicts the approximate relation (|20p . The symbols indicate results 
from the numerical simulation of 4 x 10 runs of the Langevin equation which is equivalent to the 
backward equation ((2} with absorbing boundary at 6^ bs = 0.8-7T and a starting point at 8' — 0.057T, 
(/>' — 0. The values of the other parameters are A = 0.1, p = 1, u = 10 ui r , H — and h e g = —0.1. T 
is displayed in the inset as a function of H, where H + h e s = —0.1, a = 5, ui = 10 u r , A = 0.1 and 
p = 1. Only small systematic deviations from the theoretical prediction (solid line) are visible. 

MFPT for the magnetic moment m(i) in the state fl = —p and decreases it for m(t) in the 
state Vl — p, where p = ±1. Physically, this difference in MFPTs follows from the natural 
counter-clockwise precession (if looked from above) of the magnetic moment. Therefore, for 
the two directions of the magnetic field rotation, p = —1 and p = +1, the forced dynamics 
of m(t) in the up and down states is different. Hence, the rotating magnetic field breaks the 
degeneracy between the up and down orientations of the magnetic moment for H = 0. 

Magnetic relaxation. - To illustrate the role of the effective magnetic field h e s, we con- 
sider the relaxation of magnetization in a system composed of ferromagnetic nanoparticles 
whose easy axes are perpendicular to the plane of field rotation and a > 1. The reduced 
magnetization of this system can be defined as p(t) = [N + i(t) — N-x(t)]/N, where N(^> 1) 
and iVn(t) denote the total number of nanoparticles and those that are in the state fi, respec- 
tively. Using N-i(t)+N +1 (t) = N, we obtain (i(t) = 2N +1 (t)/N. On the other hand, the time 
dependence of JVn(t) is governed by the kinetic equation Nn(t) = N^£i(t)w^n — Nu(t)wn, 
where wq denotes the transition rate from Q to — f2. This gives for the magnetization the 
well-known equation 

AC*) = + w +i) - w +i + w -i- ( 21 ) 

Because the mean residence time in the state ft equals Tq, the transition rate wq is given 
by Wq = 1/Tq. In this case, solving eq. (|2ip with the initial condition /i(0) = 1 and assuming 
H = 0, we obtain the relaxation law 

fi(t) = (1 - /Xoo) exp[-t/T] +Moo, (22) 

where r = \/(w-\ + w + i) = tq / coah(2ah c g) and fXoo = — w±x)/(w-i + w+i) = 

tanh(2a/i e ff) are the relaxation time and steady-state magnetization, respectively. Thus, the 
rotating magnetic field decreases the relaxation time and magnetizes the nanoparticle system. 
These results are not evident because they arise from the difference between the up and down 
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dynamical states of the magnetic moments. It is important to note in this context that even 
for small values of the effective field the magnetization effect can be sizable. In particular, if 
h = 0.1 and uj = 10 then \h e s\ = 10~ 3 , and for a = 50 we get \/J,oo\ = 0.1. 

Conclusion. - We have shown that a magnetic field rapidly rotating in the plane per- 
pendicular to the easy axis of a nanoparticle, lifts the degeneracy between the up and down 
orientations of the nanoparticle magnetic moment. This lifting is characterized by the ef- 
fective magnetic field acting along the easy axis in a direction that is uniquely defined by 
the direction of the magnetic field rotation. The effective field changes the MFPTs for the 
nanoparticle's magnetic moment and, as a consequence, changes the relaxation law and yields 
a net magnetization for a system of ferromagnetic nanoparticles. 
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